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^ \ In the target rest frame and at high energies, Drell-Yan (DY) dilepton production looks 

^ ■ like bremsstrahlung of massive photons, rather than parton annihilation. The projectile 

. quark is decomposed into a series of Fock states. Configurations with fixed transverse 

I separations in impact parameter space are interaction eigenstates for proton-proton 

' (pp) scattering. The DY cross section can then be expressed in terms of the same color 

■ dipole cross section as DIS. We compare calculations in this dipole approach with E772 

i data and with next-to-leading order parton model calculations. This approach is espe- 

"p^ I cially suitable to describe nuclear effects, since it allows one to apply Glauber multiple 
scattering theory. We go beyond the Glauber eikonal approximation by taking into ac- 

^ . count transitions between states, which would be eigenstates for a proton target. We 

Qj I calculate nuclear shadowing at large Feynman- jji;' for DY in proton- nucleus collisions 

>-Ch ■ and compare to E772 data. Nuclear effects on the transverse momentum distribution 

J> . are also investigated. 
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1 DY dilepton production in pp scattering 



Although cross sections are Lorentz invariant, the partonic interpretation of the microscopic 
process depends on the reference frame. As pointed out in in the target rest frame, 
Drell-Yan (DY) dilepton production should be treated as bremsstrahlung, rather than parton 
annihilation (see also 0). The space-time picture of the DY process in the target rest frame 
is illustrated in Fig. |l]. A quark (or an antiquark) from the projectile hadron radiates a 
virtual photon on impact on the target. 




Figure 1: A quark (or an antiquark) inside the projectile hadron scatters off the target color 
field and radiates a massive photon. The subsequent decay of the 7* into the lepton pair is 
not shown. 

A salient feature of the rest frame picture of DY dilepton production is that at high 
energies and in impact parameter space the DY cross section can be formulated in terms 
of the same dipole cross section as low-x^j DIS. The cross section for radiation of a virtual 
photon from a quark after scattering on a proton, can be written in factorized light-cone 
form 



d\na 




similar to the case of DIS. Here, aqq is the cross section for scattering a gg-dipole off 
a proton which depends on the qq separation ap, where p is the photon-quark transverse 
separation and a is the fraction of the light-cone momentum of the initial quark taken away 
by the photon. For shortness, we do not explicitly write out the energy dependence of agg. 
We use the standard notation for the kinematical variables, Xi — X2 = xp, r = M'^/s = X1X2, 
where xp is the Feynman variable, s is the center of mass energy squared of the colliding 
protons and M is the dilepton mass. In (|I]) T stands for transverse and L for longitudinal 
photons. 

The physical interpretation of (|l]) is similar to the DIS case. The projectile quark is 
expanded in the interaction eigenstates. We keep here only the first eigenstate, 

|g) = V^|gfea,e) + <^|g7*) + ---, (2) 

where Z2 is the wavefunction renormalization constant for fermions. In order to produce a 
new state the interaction must distinguish between the two Fock states, i.e. they have to 
interact differently. Since only the quarks interact in both Fock components the difference 
arises from their relative displacement in the transverse plane. If p is the transverse sepa- 
ration between the quark and the photon, the 7*g fluctuation has a center of gravity in the 
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transverse plane which coincides with the impact parameter of the parent quark. The trans- 
verse separation between the photon and the center of gravity is (1 — a)p and the distance 
between the quark and the center of gravity is correspondingly ap. Therefore, the argument 
of aqq is a p. More discussion can be found in 

The transverse momentum distribution of DY pairs can also be expressed in terms of the 
dipole cross section [^. The differential cross section is given by the Fourier integral 



after integrating this expression over the transverse momentum of the photon, one obvi- 



The LC wavefunctions can be calculated in perturbation theory and are well known 0, |^. 
The dipole cross section on the other hand is largely unknown. Only at small distances p it 
can be expressed in terms of the gluon density. However, several successful parameterizations 
exist in the literature, describing the entire function aqg{x,p), without explicitly taking into 
account the QCD evolution of the gluon density. We use the parameterization by Golec- 
Biernat and Wiisthoff for our calculations. Fig. |^. This parameterization vanishes oc p^ 
at small distances, as implied by color transparency |Q and levels off exponentially at large 
separations. 

In Fig. ^, we compare to E772 data [0 on low X2 DY dilepton production |^ . Most of the 
data are quite well described without any K-factor, which does not appear in this approach 
since higher order corrections are supposed to be parameterized in aqq{p). Moreover, the 
calculation in the dipole approach agrees with the next-to-leading order (NLO) parton model 
calculation at low X2. Note that the dipole approach is valid only at low X2 0. At large 
X2, this approach is not applicable and differs strongly from the parton model calculation. 
The disagreement between the data and both of the calculations in some points is probably 
due to systematic errors in the measured cross section. Preliminary E866 data |^ agree well 
with the NLO parton model calculation. 
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ously recovers ([l|). 
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Figure 2: The points represent the measured DY cross section in pD scattering from 
E772 The solid curve is calculated in the dipole approach, while the NLO parton model 
calculation (using CTEQ5M parton distributions ^TdjJ) is shown as dashed curve. The dipole 
approach is valid only at small ■ 
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2 Proton-nucleus (pA) scattering 



The main advantage of the dipole approach is its easy generahzation to nuclear targets. Fur- 
thermore, it also includes some higher twist effects that are important in multiple scattering, 
and it provides insight into the physical mechanisms underlying nuclear effects, which are 
not easily accessible in the parton model |Tl[ . 



Shadowing in DY is an interference phenomenon due to multiple scattering of the pro- 
jectile quark inside the nucleus. In the target rest frame, where DY dilepton production is 
bremsstrahlung of massive photons, shadowing is the Landau-Pomeranchuk-Migdal (LPM) 
effect. These interferences occur (Fig. |ip, because photons radiated at different longitudinal 
coordinates Zi and Z2 are not independent of each other. Thus, the amplitudes have to 
be added coherently. Destructive interferences can occur only if the longitudinal distance 
Z2 — Zi is smaller than the so called coherence length Ic, which is the time needed to dis- 
tinguish between a quark and a quark with a 7* nearby. It is given by the uncertainty 
relation, 

I =J_ = ^ (1-^)^^ (4) 

' AP- mNX2 ql + {l- a)M^ + a^m^g ' ^ ^ 

Here, AP~ is the light-cone energy denominator for the transition q — >■ 97* and qj- is the 
relative transverse momentum of the 7*g Fock state. For zi — Z2 > Ic, the radiations are 
independent of each other. 

An immediate consequence of this is that Ic has to be larger than the mean distance 
between two scattering centers in the nucleus (~ 2 fm in the nuclear rest frame). Otherwise, 
the projectile quark could not scatter twice within the coherence length and no shadowing 
would be observed. 

We develop a Green function technique 0, which allows one to resum all multiple scat- 
tering terms, similar to Glauber theory, and in addition treats the coherence length exactly. 



The formalism is equivalent to the one proposed in [12] for the LPM effect in QED. Our 



general expression for the nuclear DY cross section reads 

dlna dlna 2 J J_^ J-^^ J J 

-1 * 

X ^'^.g(a,p2) PAib,Z2)agq{ap2)G{p2,Z2\pi,zi) 

X pA(^2l)crgg-(«Pl)^7*g(«,Pl)- (5) 

The first term is just A times the single scattering cross section, where A is the nuclear 
mass number. The second term is the shadowing correction. The impact parameter is b 
and the nuclear density is pA- The Green function G describes, how the bremsstrahlung- 
amplitude at zi interferes with the amplitude at Z2- 

To make the meaning of Eq. (|^) more clear, let us first consider a limiting case for G. 
In the simplest case, the coherence length, Eq. ^, is infinitely long and only the double 
scattering term is taken into account. Then G (p2, ^2 I Pi, ^i) = ^^^\pi — P2) and one of the 
p integrations can be performed. The 5-function means that at very high energy (infinite 



5 



coherence length) the transverse size of the 7*g Fock-state does not vary during propagation 
through the nucleus, it is frozen due to Lorentz time dilatation. Furthermore, partonic 
configurations with fixed transverse separations in impact parameter space were identified 
a long time ago 0] in QCD as interaction eigenstates. This is the reason, why we work 
in coordinate space. Namely, in coordinate space, all multiple scattering terms can be 
resummed and in the limit of infinite Ic one obtains 

^frozen ^^^^ | ^ ^(2)^^^ _ f ^,,iy) ^^^^^^^ ^ _ (g) 



The frozen approximation is identical to eikonalization of the dipole cross section in Eq. (|T]). 
Thus, the impact parameter representation allows a very simple generalization from a proton 
to a nuclear target, provided the coherence length is infinitely long. 

At Fermilab fixed-target energies {^/s = 38.8 GeV for E772), this last condition is not 
fulfilled and one has to take a finite Ic into account. The problem is however, that Ic, 
Eq. (^, depends on the relative transverse momentum qj- of the 7*g-fluctuation which is 
the conjugate variable to the size p of this Fock-state and therefore completely undefined 
in p-representation. The quantum mechanically correct way to treat the in Eq. (|^) is 
to represent it by a two-dimensional Laplacian in p-space. The Green function which 
contains the correct, finite coherence length and resums all multiple scattering terms fulfills 
a two-dimensional Schrodinger equation with an imaginary potential, 



'. d At (P2) - r/^ i „ , , 

+ TTET-T, IT + nPA [b, Z2) {ap2 



G (P2, Z2 I Pi, Zi] 



dz2 2Eqa{l-a) T 

= i<5(z2-^i)5^')(p2-pi), (7) 

where 77^ = (1 — a)M^ + o?m^^. For details of the derivation, we refer to [^. 

The imaginary potential accounts for all higher order scattering terms. The Laplacian 
implies that the Green function is no longer proportional to a 5-function. This means the 
size of the 7*g fluctuation is no longer constant during propagation through the nucleus. One 
can say that an eigenstate of size pi evolves to an eigenstate of size P2 7^ Pi, so transitions 
between eigenstates occur. 

Calculations with Eqs. (|^) and (|^) are compared to E772 data [jl3| in Fig. ^ Note that 



the coherence length l^. at E772 energy becomes smaller than the nuclear radius. Shadowing 
vanishes as X2 approaches 0.1, because the coherence length becomes smaller than the mean 
internucleon separation. It is therefore important to have a correct description of a finite Ic 
in this energy range. The eikonal (frozen) approximation, Eq. (|^), does not reproduce the 
vanishing shadowing toward X2 — ^ 0.1. The curves in Fig. ^ are somewhat different from the 
ones in because we used a different parameterization of the dipole cross section. Note 
that for heavy nuclei, energy loss [0 leads to an additional suppression of the DY cross 
section. 

Nuclear effects on the gT^-differential cross section calculated at RHIC energy are shown 
in Fig. See for details of the calculation. The differential cross section is suppressed 
at small transverse momentum q^ of the dilepton, where large values of p dominate. This 
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Figure 3: Comparison between calculations in the Green function technique (solid curve) 
and E772 data at center of mass energy ^/s = 38.8 GeV. for shadowing in DY. The dashed 
curve shows the eikonal (frozen) approximation, which is not valid at this energy, any more. 



suppression vanishes at intermediate ~ 2 GeV. The Cronin enhancement that one could 
expect in this intermediate qt region is suppressed due to gluon shadowing |[T6[| . 

A nuclear target provides a larger momentum transfer than a proton target and harder 
fluctuations are freed, which leads to nuclear broadening. Note, that not the entire sup- 
pression at low qT is due to shadowing. Some of the dileptons missing at low qx reappear 
at intermediate transverse momentum. At very large transverse momentum nuclear effects 
vanish. 



3 Summary 

We express the DY cross section in terms of the cross section for scattering a qq dipole off a 
proton. This is the same dipole cross section that appears in DIS. At low X2 and for proton- 
proton scattering, calculations in the dipole approach agree with calculations in the NLO 
parton model. Some E772 data points are not well described by either of the approaches, 
which is probably due to a systematic error in the measured cross section. 

At very high energy, the dipole approach is easily extended to nuclear targets by eikon- 
alization. At lower fixed target energies (E772) the eikonal approximation is no longer valid, 
because the size of a Fock state varies during propagation through the nucleus. Therefore, 
transitions between interaction eigenstates (i.e. partonic configurations with fixed transverse 
separations) occur. 

We develop a Green function technique, which takes variations of the transverse size into 
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Figure 4: Nuclear effects on the DY transverse momentum distribution at RHIC and LHC 
for dilepton mass M = 4.5 GeV and Feynman xp = 0.5. 

account and resums all multiple scattering terms as well. For light nuclei, calculations with 
the Green function technique are in good agreement with DY shadowing data from E772. 
For heavier nuclei, also energy loss becomes important. 

We have also calculated nuclear effects in the transverse momentum distribution of DY 
pairs at RHIC energy. The DY cross section is suppressed at low transverse momentum. 
The expected Cronin enhancement at intermediate qt ^ 2 GeV is reduced because of gluon 
shadowing. Nuclear effects vanish at very large qt- 
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